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Rung<fKutta Mothoda, Hw'.apioal Analyoio, Opdinary Diffepontial 
Squationo 

icalod Runao-K utto Algorithms for HomUlnq Donoo Outftut 
Summary 

Now, low ordor Runga-Kutta algorithms aro dovolopod which dotor- 
roina the solution of a system. of ordinary differential equations . 
at any peint within a given integration stop, as. well as at the 
end of each step.-4?he new, sealed Rurtgo-Kutta methods are designed 
to be-used with existing Rurtge-Kutta formulas., using the deriva- 
tive evaluationa of these-definlng algorithms. as the core of the 
new dystem« Thus, for only a slight increase in computing time, 
the solution may be generated within the integration step, improv- 
ing the efficiency Of tho-Runge-Kutta algorithms, .since the step 
length need no longer be Severely reduced to coinSide With- the de- 
sired output point. Scaled Runge-Kutta algorithms are presented — 
for orders 3 through 5, along with accuracy Comparisons between 
the defining algorithms a nd the lr-scaled versions— for a test pro- 
blem. 


Kunge^Kutta^Vepfahpent Humepieaha Mathematik, GewShnliahe Diffe- 
ifentiatgteiohungen 

Skalierte Runge-Kutta -Forme In fUr dlchte LOsuncsausaabe 
Obersicht 


Neue skalierte Runge-Kutta-Pormeln dritter bis fttnfter Ordnung, 
die die bttsuttg eines Systems von Dlf ferentl algleichxmgen an be- 
stimmten-Punkten Innerhalb und am Ends eines jeden Integrations- 
sohrittes bestimmen kiJnnen, werden hergeleitst. Die neuen skalier- 
ten Algorithmen werden unter Verwendung der klassisohen Runge- 
Kutta-Fomteln entwlckelt und bentttzen deren Schema von Funkbions- 
auswertungeh. Auf dlese Weise kann die IiOsung innerhalb des Sohrit- 
tes roit wenig mehi Rechenzelt bestimmt werden, so daS die Lei- 
Btungsfdhigkeit der Runge-Kutta-Methoden in hohem Grade verbessert 
wlrd. Die Qenauigkeit der Lttsungen wird mit derjenlgen der klassi- 
schen Forroeln an Hand eines Beispieles verglichen. 
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1 . Introductio n 


The Rungo'Kutta (RK) algorithm, dooignod to tsdat tho initial val- 
ue problem 

(1.1) “ f(t,y) , y(to) ° yo 

whoro y is an n-vootor, is an explicit, singlo-stop algorithm 
which uses a linear combination of derivative evaluations to ap- 
proximate a truncated T^/lor series. The high accuracy and the 
simplicity of the algorithm lendo itself readily to the solution 
of many classes of ordinary differential equations (ODEs) . To per- 
form efficiently, however an RK method must operate with an opti- 
mal step size, consistent with the specified truncation error cri- 
teria. If frequent data output is requested, e.g., in iterative 
schemes to determine the perigee or apogee of an orbit, this ef- 
ficiency may be reduced significantly. This difficulty arises be- 
cause the Runge-Kutta algorithm generates the solution of tho ODE 
only at the end point of the integration step,_thus requiring step 
size reductions to output-data at specified values of the indepen- 
dent variable. Shampine, Watts, and Davenport [4], however, have 
shown that such adjustments to the step size seriously impair the 
efficiency of the method. The Runge-Kutta algorithm must operate 
using near optimal step sizes if the method is to be competitive 
with other numerical methods for solving ODEs, and, thus,^jan al- 
ternate technique for treating the problem of frequent data output 
is essential* 

Scaled, Runge-Kutta algorithms may be developed which determine 
the solution of the ODE at any point t* = t* +o h,. once tho solu- 
tion has been evaluated at ti = to +h, where the scaling factor, 
o, generally ranges between 0 and 1 . By using the derivative eva- 
luations generated to obtain y(ti), this now scaled solution may 
be determined at a cost of, at most, only a few additional deriva- 
tive evaluations. Thus, for a given RK algorithm, a family of 
scaled algorithms may be developed, which, when used in conjunc- 
tion with this defining algorithm, gives tho solution at any point 
within tho interval of integration. 
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PAtSS 

bwnk not Wr.MRn 


The tlorlvatlon of tho ooaloa algorithm lo givon for oraoro throo 
through five, boirtg applied to oooffioionte of defining algo- 
rithran dovolpod by G. Fehlberg (1] and D.Q. Oottio (privato eom- 
munloation) • Nutnorloal rosulto« comparing tho accuracy of tho 
Denied oolutlon to the accuracy of tho dofinlng algorithm, oro 
prosontod. 


2. The Scalod Runqo-Kutta Algorithm 

Tho Rungo-Kutta algorithms of order p, used to solve Eqn. (1.1), 
assumes the form 

r 

(2.1 *a) y « yo +h Z C. f . 

j«0 J J 

whore 

(2«1«b) fo ~ f (to # yo ) 
and 

(2.1.0) f . = f(to +a.h, yo +h 

for j «1^2#..*,r^ where the a, 3# and C coefficients have been se- 
lected so that the algorithmic solution# y# is equivalent to a Tay- 
lor sum of order p* Specifically# a term by term comparison of the 
Taylor series expansion of the algorithmic solution to tho Taylor 
series expansion of the solution vector results in certain trunca- 
tion error coefficients# TEC^^ ^ #» i.e. # expressions in a# 3# and C 
for each order# i# with j ®1#2#..*#^^ where increases with in- 
creasing order. For an algorithm to be of order p all TEC^ . must 
vanish for i«1#2#...#p) j » 1 #2# . ♦ . #Xj^* Those vanishing trunca- 
tion error coefficients are referred to as equations of condition. 
(Soo Table 1.) 

&j #0 determined from tho relations 
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J"1 

); 

k o 




J.k 


■■ « 




for J 1,,?,. 


. . ,r. 


Till.' lorin or Uk' iioalod ulyorilhm parallolo that of tho Uofinlny 
al.jorlLlim olnoo tho ofCicluncy of ttu> now mothodo roijulroo tins 
uuu of u.o chrrlvativo evaluakiono (2.1.b) and (2.1.c) an tho eoro 
of U.o now avutom. onu auaumus that tho aolution has boon advanced 
from i„ to ti '=to +h, l.o., that fi, fac..,fj. have boon dotormin- 
od, and that tho aolution la ro.juirud at oomo Intcrmodiato point, 
t* to 1 11 * -^to +0 h, with o a scaling factor, where generally o e 
(0,1). More uiJoolf Ically , 

(2.3) h* -oh 

Tho scaled algorithm assumes the basic form of the defining algo- 
r i thm 

r* 

(2. 4. a) y(to +h*) ^ yf ~ yo -fh* S Ct fj 

j«0 ^ J 

where* 


(2*4. b) « f(to, Vo) 

and 


. 1-1 


(2.4*c) f* St f(to +a* h^i yo +h* S 3? u f(!!) 

k«0 ^ 

for j -1,2,...,r*, whore r* 5r. If a* and 3* are chosen so that 
(2.5) fj-f^ 

for J -1,2,...,r, then the major expense of applying tho algorithm 
has been avoided since fi, fa,...,fj, have already boon determined 

during the evaluation of y(t,) . Thus, only f^ f need to bo 

evaluated. Comparing eqn. (2.1.b) and (2.1.c) with (2.4.b) and 
(2.4*c), iron|;>oct ivoly, tjivos ~ f ^ if 
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(2.c,n) a (tjh/h* r. a^/a 
and 


(2.6.b) o 0^^^h/h* a 1^/0 

for j “'Ij2>...,r) k ^0(1,2»...)j“1. 'i'hui.i( many or iIk‘ cocri loJi.'it( u 
of tho ocalotl alijorithni aru produtormlnoU, bu.lu<i i-claii'U i<i iliolj- 
corroopondlng cocfficiontw in tho dufiiuj uliiorltlim by a (aotor of 

1 / 0 . 


Tho solution of the sealed algorithm dopundu upon Uu> aimumpt: ions 
imposed during tho derivation of the defining algorithm, I.e., 
those assumptions must generally hold in the sealed system in 
order for the new system to satisfy the same equations of condi- 
tion. Tho most general assumptions imposed during tho d<’Vclopment 
of the low order RK methods concern tho terms 


(2.7) 


■j»i 


j-1 

£ 

k=1 


^j,k®k 




sincG setting 


( 2 * 8 ) 

for j causes a number of equations of condi-* 

tion to become redundant* Clearly, after substituting tho scaled 
parameters, Oqn* (2*8) is also valid in the scaled system* An im** 
portant consideration is the analysis of the higher order trunca-* 
tion error coefficients and tho stability properties of the algo- 
rithm which are tho major criteria used in the selection of the 
coefficients for the defining algorithm* In general, during the 
solution of the equations of condition, several a, 3, or C para- 
meters remain as free parameters which may be t^olected to reduce 
the TEC of tho next higher order improving the accuracy of tho 
method* Theses same parameters govern the stability properties of 
the method, i*e*, they determine the error amplification or damp- 
ing due to tho algorithm* Eof the linearised equation, y* ~A,y, the 
error amplification factor, r(Xh) , is determined by the polynomial 
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(Xh)'- , 


III'- I 

' n v.(Ali)' 

K) k*orclpr+1 


will !«• Iim M|, imiiw In «, n, mid c, and wliorn li in 

**''' '**"!’ '■Ic'!'' iitid A 111 mi I' I iii'iival IK' (,if lIk.i i)iu.'ol.ilan iiiali'lx« 
('l/iyl, I'V.i I ii.'it I'd ill I ■ l.o , 'i(Ali)|i 1 di'fliu'ti Uk.i ij|;al.».l 1 ity bound* 
my tor ii .jjv.ii Id! aliiorll Inn, Ki'itu'iiiboi'inij that li will Ixi ooalod by 
< 1 , III.' Till' ii riiir. mid Um iitablllly ijolynoiiilal will altio bo affoot- 
' d niijni r Jomii ly by thin ncallni) paramo* or. 


llii' cxpoiino oi ii|,i lyiim tlio ucalod alttorlthm In dotorminod by tho 
luiiiibcf <j( additional doi-lvatlvo ovaluatlonu roqulrod. {'or tho third 
ordor mothud iiLudlod. no additional dorlvatlvo ovaluatlonn aro 
iioodi-d, I.O., tlu! uolutlon may bo obtained at any intormodlato 
valai- of tho Indopondont variable at a coot of only a few arlth- 
iiiotio oporationn uood in generating the C* coef flclonte. Using one 
uddli tonal derivative evaluation, a fourth ordor solution may bo 
determined anywhere within the integration step. After this addi- 
tional evaluation has been generated, the solution may be computed 
at furtlKV values of t at a coat of a few arithmetic operations, 
the co.npiexlty of the solution of tho fifth order, scaled algorithm 
dictate!) the necessity of two for,as of such algorithms. For two 
additional derivative evaluations, a fifth ordor solution may bo 
det)'rmlnod at any point within tho integration stop (Horn, t2)). 
Kach further solution, however, requires two additional derivative 
evaluations. Thus, for dense output within a given step, a second 
form of the scaled, fifth ordor solution is needed. Using five ad- 
ditional derivative evaluations, a fifth ordor solution may be ge- 
nerated within the Integration step with the added advantage that 
further solutions may be determined with no additional derivative 
evaluation expense. Such a formulation is of particular importance 
In iterative schemes involving an RK integrator. 
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3 1 A Pumin- Huttn ,fyii[orj tiim 

UBtny fl third prUnr, Rumfu-ltutta (ildpritiun (due u« P.n. hemn, 

0 «p Tahlo 2 id)) having an t'mljcddnd tincond ordiT iiulut Iuh fur ntc|. 
oi«« eentrel, a (luuU'd, lldrd orUor a]i|uritlim r y i.o d<'Vi'Io|'<'d 
having eemparativoly iow truneatlun im'ur ecK>. . UMc niu throuiihtiut 
feho interval nu well au reanonahle utahUity i'r.i|..i-rUi-n. lilnee the 
defining algerithro^ t^K'l'(2)3| Itao the advantage that the final tie* 
rivative evaluation for one utep egualu the InUial evaluation for 
the oubooquont otup, third order uolutionu may be evaluated at any 
value of the indupendont variable at a ooiit of tliree evaluatlenu 
pur integration utup (pluu the orithmotic operatlrmo rt^iuired for 
thu C* coeff iciuntu) f unlcuu the uolution it) to be uUvanoeU from 
0 noalud uolution evaluatlonf in which caue one a^d ional • rl- 
vativu evaluation is needed. 

The equationu of condition for the third order uiale-’ j ..m 

are 


(3.1) C8 +Cf aij +CJ aj tCf oV(j + 1) 

for j “ 0 , 1 , 2 , where cto ^0 and 


(3.2) Cj 0a,iOi +CS(3a,,a, +33,aOa) ■= o“/ 


Whore the and 3^^,^ coefficients are from the defining algorithm 
(Table 2 . a.), while the C* ’coefficients arc the weighting factors 
in the scaled algorithm. These same equations detormlno the coef- 
ficients of the defining algorithm, except that O 3 -O and the 3 
coefficients are defined by ^ 3 ^)^ where the Cj^s may bo written 
solely in terms of the o coefficients. The C* coeff lei onts of the 
scaled algorithm may then bo written! 


(3.3) 


Cf 


= “0(2o -3) (2 -3aa) 
6ai (ai -oa) 


- 0(20 - 3)/3 


(J..1 Cl . rw q _ -,l(2 

6aa(aa -ai) 

and 
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(3*5) C% c=ci»' a 

v.»l th 

(3.0) 1 - (C< +CJ +CS) ° (9 +0(-12 +5o))/9 . 

Table 2.b presents values of GS, the Euclidean norm of tbo fourth 
order truncation error coefficients, 

(3.7) GJ n.. ( E T* =)’/2 ^ 

J^cr 1 ^ 

whore ti.e order of the TEC terms i =4 and where the number Of 
those terms for order i,X^ = 4, with the Tj,^j^ coefficieftts- (given 
in Table 1) , as well as tne stability limits on the real and im* 
aginary axes, R»- and I*, for various, values of o. Results compaes 
ing the accuracy of the scaled algorithm to that cf the defining 
algorithm for the problem in Section 7 may be found in Table 2.c. 


^ • A Scaled, Fourth Order Runge-Kutta Algorithm 

The fourth order, Runge-Kutta algorithm due to E. Pehlbcjfg (1J^ 
having an embedded fifth order method for step size control, has 
particularly desirable properties for scaled algorithms. More-spo 
flcally, once one additional derivative evaluation has. been made, 
the solution may be determined at any number of points within the 
given step for a few arithmetic operations, while still maintain- 
ing the accuracy of tlie defining algorithm. 


Since eqn.(2.8)ls imposed for i =1 and 2, the only remaining equa- 
tions of condition reduce tot 


(4.1) S C* = oV(j+1) 

If trn ^ 


k*^0 


for J “0,1, 2^ 3, where ao «0, Cf =0, and 


(4.;!j E 1 “ O 

k«=2 ^ ^ 
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far j ■0, SUffiolant puramotota exist to solve Eqn, (4«2) for 
j a fifth order term, .which roaults. In overall lower TECs. 

Since several solutions are sought within the stop, the M , co- 
offlclent may not be used to solve Eqn. (4. .2).. Instead, 6nly the 
C* coefficients are available for generating the scaled solution 
with the for o«1, determinod form Eqn. (2.4). and the ro- 
sining related to by a factor of l/o for k =0, 1 , 

(At first glance.,, a solution appears possible for r *5. Eqn. (4. .2) 
however, would require infinite c* coefficients.) * 


fiqn. (2.3), for j «1 and. 2, determines the scaled p« and. p*' 
eoefficients with the unsealed ^6 3 free pa^- 

roetors, givino m- o,4 6,3 

3 


(4.3) 


^t.± 


a|(2a^-3a.).6E c^(a^.-a,) 


k°1 
60 a» (d4 


• 081 ' 


with p |^5 determined, by permuting Indices 4 and-5, where ♦ denotes 
coefficients from the scaled algorithm, while the absence of ♦ dc^ 
notes coefficients from the defining algorithm or a* or p* for 
0 -1. Eqh. (4,.1),. for j =1,^4., determines Cl, CS, and Ct!in terms 

of Cl, ,c«, and o-.feqn. (4.2), fOr j -1 and 2, then gives eland C| 
solely in terms of a. 


(4U.d) Cj - (AAa..±DDaCj +EEaCj) 


Where 

3 

(4.4.b) AAa - S Bp- nO** •> 

Ii-l — 

= (03/4 -oa(a» +d4)/3+c oa 04 Jvtcia (o# -a») (a* -o,)) 

(4.4.C) DD*. . -o«(ai -as(a3+d8)..*aao*)/(aa(aa -oj) (o* -o,J) 

(«.«.d> EEa • -06 (a| -a8(da+d,) +Oafll8j/(tta (a, -o,) (a* -a*)) 

with CJ ahd CS determined by permuting indices, 2,3,A...ithe..C| and 
• tarma become 
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i 1 . «- bi8<l3 . . 

('I . «»«\) L’#-, - J .J1 qJ 

j Di lD,»;< ** b.»^ b , 2 

. ^ ) >1 1^2 i - b;. i<t 1 j j 

(■I. .!■- (•,’ >: r-K-* •-, -.--io 

j; ^ bl lb;>;. - b;. lbl 2 

win *l 


b, 


i-1 


.1 ' Ji2 «k \,i 


' - .«. IV 


l.=x h -k ’ '\-,i 


< .i t J( t i 


•= J2 Pk.i 




Wh. re DUn '^1, DDo =0, lit:,. =0, and EEe = 1 . Then all C* coefficients 

tiiv.f.‘» jn' Vf) .lt,iun aa polynun;.ials in a* 

Unjny the RKF(/])5 coefficients listed in Table 3*^ and chosing 
Ufo a^if and Ps/ i “Pe/a-bou ^ O, the scaled coefficients become 

(4,6*a) Co - 1 " o (301/120 + a (-269/108 +0(31 1/360) ) ) 


(4.0*b) •'* ^ 0(7168/1425 +o(-4096/513 +0(14848/4275))) 


(4.<»*c) cl u( ;’8bi>1/8360 +0(199927/22572 +o(-371293/75240) ) ) 


(4.b*il) ci - 0(57/50 +o(-3 +0(42/25) ) ) 


(4*6*4.»J i** o{-96/55 4 o(4u/11 +o(-l02/55) ) ) 


(4.6.i) -M 3/2 +o(‘ 1 4 0(5/2)) ) 


15 


with C? 


0, and. 


(4«7-»a) 06/6 » 31/6 *® 1/(6 o) 

ahd 

(4*7J;^ 36/.S — s— i/(3 o) * 

As a approaches unity, all fifth order TEC approaoh zero,^givind a 
true, fifth order solution at o »1. one should note that since the 
P coefficients are multiplied by o during the evaluation 

Of fe, fj remains unchanged When differing values of o are used. 
Therefore, fs needs to be evaluated only once.. 

Coefficients for the defining algorithm, RKP(4)5 are listed in 
T^e 3, With the values of cf, Eqn.(3.7), the Euclidean norm of 
^e fifth ordef truncation error coefficients, as well as the sta- 
bility limits along the real and imaginary axes_being presented 
In t^le 4. a for various values of o. Accuracy comparisons between 
the fourth order solution and the fourth order scaled solution are 
given in Table 4. b^ with results for the fifth order solution hav- 
ing a scaled, fourth otder solution being given in Table 4.c for 
the problem stated in Section 7 , 


?ive-Stage , Fifth Order. Scaled RK AloQr^^^h«. 

Due to the complexity of the equations of condition for the fifth 
order, Runge-Kutta method, the generation of a fifth order, scaled 
algoritto requires more additional derivative evaluations than are 
needed for the lower order, scaled solutions. To generate a scaled 
solution Which is valid for all values of o. once the set of addi- 
tional evaluations has been made, requires five more function eva- 
luations (five stages) . Such an algorithm is Ideally suited for 

the problem of dense output, e.g., for Iterative schemes over a 
given step. 
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The equations of condition for the fifth order RK algorithm are 
given in Table 1, with additional assumptions, Eqn, (2.8), being 
used to force some of these equations to become redundant, fiqn. 
(2.8), for j >=1 and 2, reduces the equations of condition to 


(5.1) 


j 




aj = oV(i +1) 


for i =0,1, 2, 3, and 4, 
r* 


(5.2) -(5.4) 2 CJ ai 3 . . = o 

joi J J 3»K 


with i =0 or 1 for Ic =1, and with i =0 for k =2, 


(5.5) 


r* j-i 

1=2 J k=1 ^ 


o‘*/20 


and 


(5.6) 


r* j -1 

S 2 B 


j,k r.o 


j=4 J k=3 

where r* =10 and where ♦ indicates coefficients in the scaled al- 
gorithm, while the absence of * Indicates either coefficients in 
the defining algorithm or a^, PJ for o =1. 


The solution of (5.1) requires the inversion of an nxn Vandermonde 
matrix (i.e., a matrix with element ® con- 
stant) , whose analytic solution involves the coefficients of the 


Lagrange polynomial, with the coefficients given by 


(5.7) L^(x) 


ft* 

k=1 


(X -u^) 


k?i (Wi -uj,) 


n 

Z 


j»1 


^i.j ^ 


j-1 


on the defining set M=(ui , . . . ,u_) where ' indicates k -"i 


(i»j) element of the inverse matrix. 


and wher« 


A scaled, fifth order algorithm in which the solution may be de- 
termined at any number of points within a given integration step 
at a total cost of five derivative evaluations (plus the arithmetic 
operations used In evaluating the new C*-*s for each solution point), 
may bo derived by setting the C., coefficients of the existing f^ 


terrta. equal to zero. Then, by adding five new derivatiAre evalua- 
tions, £J, f?, fj, and.fto, With 3} , »0, for 

16,. and iiobosing Eqn. (2.8) , .^i = 1 ,1' and 3^ the remaLing ’ 
equatiottd-of condition beOoae^ 1 )Eqn. (5.1) „ for i - 0 , . . . , 4 , 
and 2i Eqn. (5.6). Eqn. (5.1) determines all c* coeffioients’in 
tanas of one C*, say C?,. while Eqn. ( 2 ^ determines 3f ,, Sf 

^^5 5. 3 ^ 

free parameters, AS in the fourth order, scdiel algorithm all 
^j,k '^’oe^^icidnts are_related to their corresponding 0 . . for- 

0 - 1 , end-therefore, fj t*, need to be evaluated oJiJ once 

for,a given integration step. The c* coefficients may be written- 

^5+6 “ /(( 6 -k)aj^g))-C|(^fi^at"’' ®j, 5 -k/®j+ 6 ’ 

with j-1 ,2,3,4, where the-a^^,^ terms are defined on the set 
M - {o»,« 8 ,a»,o,o) , — Eqn. (5.7), and the 0^^,^ co efficient s 

(5.9) 35^,^ = ^“”/(5-m)-KK^^,.J V2,4-m/«^ 

for j-6,.7,8,9,10, k»3,4,5 where a,^. 2 ,m <3efined on the set 
« {aa^oUft/as}# fiqne(5*7), and whdr^ 

(5.10) KK, ^ 3 

for j®7g8^^#lO wltih KKe^^Oa Now defining 

«.d .olvin, (S. 6 ), CS b. vrltfn a. a polyaoMal in o, 

10 4 


(5.12) C| « 


(2 AA, .SS, 0 *^) 
i*>7 k»1 J'* J 

To" 

SSe + E BB 4 SS 4 
J-7 5 J 


from Which all eg terms may be written as fourth order polynomials 
in Of 


18 


(5.13) CJ = o(Aj j 

for 1»6,7,8,9,10. where Cj=>0, for j=>1,2,3,4,5 and whofe 
Cji“1-(C8+C^l+C8+Cj+Cij<») . 

Coefficients are listed in Table 5. a, for the defining set, 

RKJ'(4)5, (Table 3), while the nonh of the sixth order TEC terms, 
GJ, Eqn. (3.7) and the stability limits along the real and ima- 
ginary axes are listed in Table S.b. Test results comparing the 
accuracy of the scaled solution to that of the defining algorithm 
are given in Table S.c for the problem stated in Section 7. 


6. A_Tw> jtage,. Fifth Order, Scaled RK Algorithm 

The RKF(4)5 algorithm has a fifth order, scaled solution available 
for only two additional derivative evaluations, (Horn, [2]). (Each 
additional solution, however, requires two feors function evalua- 
tions.) Efficient use of such an algorithm may require a priori 
knowledge of the o values, since dynamically generating the scaled 
coefficients could use more computing time thaft the actual deriva- 
tive evaluations.. Thus, one may need to use o values for which the 
scaled coefficients have been stored. By adjusting the step length 
externally so that the desired output point coincides with a pre- 
determined value of o, the Runge-Kutta algorithm may be operated 
at near optimal step length, while still determining the solution 
at the desired output point. For dense output within a given step, 
however, the five stage algorithm _ (Section 5) , is more efficient. 

While a single stage, scaled solution appears possible from Eqns. 

(5.1) -(5.6) , such a derivation leads to infinite C* values (for o*1), 
so that one must Impose two additional derivative evaluations. Then 
Eqn. (5.1) may be solved for j«1,2,...,6, giving 

(6.1) a^,7.^o’“V((8-k)a^^,) 

defined upon the set. M=(aa,a 3 o,) for J=1,2.,...,6, where Ct«0, 
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and whore CJ Is dotcrmlnod from 


(6.i) 


cj “ i-(ct+cj+cs+c5+£j+c8+c?[) 


Eqns.. (5.2)-(5.4) determine pj,, . 0 ?,, , and 3?, a with 3S,a 
parameter, giving 


aa a froo 


( 6 . 3 . a) 38 ^ 1 = (i<i-Koa,)/(c*j 7 (tie-a,)) 
and 


(6.3.b) 38^,0 (Ki~Koa,)/(c8o(a,-aJ) 


whare 


(6«3eC) C*OL^ 


for j«0,1a and 


k «2 




6 

S 

3=2 


(6.4) 38a “ -( 2 ^8) . 


With 3J^^ as an additional free parameter, Eqns. (2.8), fe»r i=1 
and 2, determine 38^3 and 38 „ , 

( 6 . 5 ) , , (Ka-K,a„)/(o tt,(aa-a„)) 

aftd 


(6. 6, a) $8 „ . (Ka-K,a3)/(o tta(aa-aa)) 

vrtiere 


(6.6.b) 


4*1 


/ (i+1 ) -3e, , 4-3#, a oa-3«, s 4 


for i=1,2. Eqns. (5.5) and (5.6), along with Eqns. 
j"/, determine the remaining unknowns, bi n* 
Bqn. (5,5) gives '** * 


2.8 for i=1 ,2, 
B8,s and 38^«. 
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(C.7) 


and 


6 

P7»a ^ (o<*/20-^S^ CJ Pj^/C$ 


9 


(6^ <= CJ V , 

gives 

6 5 

- 2 S . - S A. p. , 

(6. 9. a) ^ 

°"j 23 °1 ®).t>*“7,ll 

whores 


(6*9. b) As a (Ka-Ka(a4+(XB)+Kia4a8)/(aa(a3-a4) (a^-as)) 

(6*9*c) Ba a -ae (<Jt|*<Xs (ct4+aB)*fa4a8)/(aa(aa-a4) (aa-a») ) 
and 

(6.9.d) = P7^i-P7»i<ii-e7.aaJ 

With A», B., As, and Bs determined by permuting indices 3,4, and 
5 on the a terms. The remaining may be determined from 

(6.10) - Aj^+B^ 

for k a 3,4,5 . 


The free parameters, aj, a?, pj^^ and , may be used to deduce 
the sixth order TECs for a given value of o. The BKF(4)5 coeffi- 
cients, Table 3,_have sixth order TECs which are quite sensitive 
to the choice of aj and a|, but, in general, only mildly sensitive 
to the choice of PJ^, and pj^^ . in fact, one may sot the free PS 
terms to zero without greatly changing the resulting norm of the 
TECs, GJ, Eqn.(3.7). Sucjgestcd values of oj and o|, and the re- 
sulting GJ are presented in Table 6. a for a wide range of o, whJ.lo 
several sets of scaled coefficients are presented in Table 6.b. 
Test results comparing the accuracy of the scaled, fifth order so- 
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lution to that Of tho defining algorithm aro givon in Tablo 6.c 
for tho tost probloro stated in tho following oootlon. 


7. Toot PrOblom 

The two-body problem with an elliptical orbit of ecoontrlcity 
0.6 is described by tho following system of eguationst 

(7.1) y{ = ya, yi o y^, y. „ .y/,.3 

whore r = (y?+y|J"». and whore y, (O) = o.4. y,(o) = o.O, y 3 (o) « 

Keow" “ gonerated^y solving 

Kepler s equation, u - 0.6*sir(u) - t. for u giving y, (t) = 

COS(U).0.6. y»(t) O0.8*sin(u). ysCt) » -sln(u)/(1-o.6*cos(u)) . 
and y a i tl-v O.0*coa(u)/(1-o.6*cos (u) ) . 

The solution Of Egns. (7.1). has been generated for several speci- 

and 

oJT K ' SKFU)5 fifth order solution. This 

^age has also been adapted into an RKT23 routine, which uses the 

Zttllr generating a third order solution. Both 

rllZTilZTZT Z"" " etep-by-step mode with the 

reacting solution for each sten being scaled for various values 

of o AD..computations were performed in double precision (16 di- 
gits) on the AMBAHI, 470 V/6 and IBM 3020 computing system. 

^ - the oev^b. Lr- 

The average of the absolute value of the errors of each component 
tZ ll <0 <t<2n). are presented in Tables 2.c. 

ordels fd solution. 4.b and 4.c for the fourth 

ution. and-6.e for tho two stage, fifth order, scaled solution 
e e e defining solution has remained unaffected by the scaled' 
olution. The solutions for the third order integral al iLtl 
for requested accuracies of l.bO-02 and 1.0D-04. whereas 11 1 
fourth and fifth order solutions aro listed for I.OD-04 and 
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1.0D-06, oinco, In gonorol,. ono would oporato tho two mothodo using 
difforont spool flod tolurancos. Tho tabloa show tnat tho soalod 
solutions maintain tho accuracy of tho doflnlftg algorlthfti. Of par- 
ticular importanco ai*o Tabloo 4.b and 4.c, which compare tho fourth 
order, acalod solution with tho defining algorithm being oporatod 
a fourth and a fifth order method, respectively, Tho RKr(4)5 un- 
sealed solution in Table 4,c (of fifth ordor) is used to generate 
a fourth order sealed solution. Thus, the fourth order errors in 
tho acalod solution result from the local trancation errors and. 
not from tho conditions at the beginning of the integration stop, 
which are of fifth ordor accuracy. Since the error estimate re- 
quires that tho difference between tho fourth and fifth order so- 
lutions be less than a prescribed tolerance, one would expect ex- 
cellent agreement between the fourth order scaled solution and tho 
fifth ordor defining solution. Thus, for determlnlne data through- 
out a step, the fourth order solution may suffice. If, however, the 
solution is to be advanced repeatedly from a scaled output point, 
the user may wish to generate a fifth ordor scaled solution, since 
frequent use of a fourth order solution may tend to degrado-tho ac- 
curacy of a fifth order ^ defining algorithm* 


8* Conclusions 

Scaled Runge-Kutta algorithms which generate the solution of an 

ODE within a given integration step have been developed for Orders 

three through five. Each scaled algorithm is used in conjunction 
with a defining Runge-Kutta algorithm with the derivative evalua- 
tions of tho defining formula forming the core of the new scaled 
methods. For the third order, scaled solution, no further deriva- 
tive evaluations are required, while for the fourth order, scaled 
solution, any number of scaled solutions may be generated within a 
given step at a total cost of only one added derivative evaluation 
(plus tho arithmetic operations involved in generating the scaled 
coefficients) . A scaled, fifth ordor RK solution is available for 
an additional two derivative evaluations, but each further solu- 
tion requires two more such evaluations. This two stage, scaled 
algoritlim may also require storage of tho RK scaled coefficients 
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parametor er « i-pUpf knewloaqo 

af tho 4o8irod output point if tho eomputlng time ier the eeaioa 
ooeCCioionto oqualn that noodod for derivative ovaluationo. A more 
gonoral fifth order, oealou algorithm pormito any number of eolu- 
tiono to bo gonoratod during a givon atop for a total of five addi- 
tional dorivativo ovaluationo, making tho method moro eui table 
for donoQ output within a givon otop. Toot rooulto ohow that tho 
accuracy of tho ocalod algorithmic oolution io comparable to tho 
accuracy achiovod uoing tho defining algorithm. 
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Tablo 1. Tnmoatlon litror Cooffioionto, ,,, fot 'um.|o-KuUa Al- 
gemitroo of ofdor i, lni,2,3,4,5. 


"m” Jo 

" 2.1 “ I “1 - I 


"3.1 “ I ( C, aj - 1, 


_ r i-l 

“ J 2 ‘ 5 

"4.1 “ 1 C, al - 1, 

" 4,2 “ ^^2 Ji “^l.j ^ 


r i-l 
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’’*•* ■ 1^3 jSj “l.J jS, °3.k "k ■ sV 

’S,! - JJ ij, =l «} - 5 ) 

■'5.2 ■ 5 =i »? ll, h.i »J - 1^) 

’'5,3 ■ i 1,2, =1 «i ‘2I h.i - A' 
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m ^ j -1 

^^5.4 “ ^ *^1 ®i ® P 4 4 S Ba , a, - 
i=3 ^ ^ j=2 k=1 ®k 


k=1 •" 

^3 _ 1 


’■s* ■ i '.i, =i ‘i’ «5 - fo' 

n, ^ j-1 

^5,6 * 2 ^1 H S 3 . * a, L 

i=3 ^ j*2 3 Jc«1 ^ 40 

m 1 ^ i~l j-1 

^ i=3 ^ j«2 k«1 ^ ^O 


r 


i-1 j-1 

k-1 

= £ 
i=4 


S 3, . S 
j=3 ^'3 k=2 

“j,k if, \,i 

-^ 

r 

i-1 

■ 55> 

£ 

i«2 

Ci ( E 3. . 
^ j=1 ^*3 
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Table ? n Coefficients for a Runqo-Kutta (Second) Third. Order 
Algorithm, RKT(2)3 


■ 


e 

j»k 

C 

j 

c 

j 


1 

1 


u 




1/2 

2/9 

■H 


1/2 



0 

1/3 



0 

3/4 


0 

4/9 

■ 

■ 

2/9 

1/13 

— tt/g 

1/2 

0 


Hotfej 3rd Oder salutiohs ate formed *itb coeCd^iciehts;. 
2nd order 8om. t i.j a o are formed kU4j coefficients. 

t 

= *f) I 

j=0 ^ ^ ^ 


Error astimato « ly • yl 




Table 2.b Characteristic Properties for the Third Orddt, Scaled 
RKT(2)3 Algorithm. 


0. 4liJ1l!)9n-01 

0. SJ9336tfD-03 
J. 1 }9»873n-i)2 
0.39587960-02 
0. •> *25054 0-02 
0. 94338575-02 
0. 12566250-01 
0. 13793420-01 
O..U02.6080-01 
0.2J187700-01 
0.2*21 84.7 0-01 
0. 11072990-01 
0.3-y71779D-01 
0.31128375-01 
0.31283870-01 
0.37173570-01 
0.31773590-01 
3.J100644 10-01 
0.41019790-01 
0. 41610610-01 


(-2.512, 0) 

(-8.056, 0) 
(-5-J681, 0) 
(-4.648, 0) 
(-4.046, 0) 
(-3.644, 0) 
(-3. 354.^ 0) 
(-3.136, 0) 
(-2. 966,, 0) 
<-2.831,_0) 
-2.722, 0) 
(-2.634, 0) 
(-2.564, 0) 
(-2.508,. 0) 
(-2.466, 0) 
(.-2.436, 0) 
(-2.418, 0) 
(-2,414, 0) 
(-2.424, 0) 
(-2.454, 0) 


to,. 1.732) 

(0* 8.247) 
(0, .7.777) 
(0, 4..«72) 
(0, 4.008) 
(0* 3.551) 
(0„ 3,211) 
(0, 2-.945) 
(0, 2.729) 
(0, 2.550) 
(0, 2.399) 
(0>. 2.270) 
(0, 2.158) 
( 0 , 2 . 062 ) 
(0, 1.978) 
(0, 1.907) 
(0,_1,847) 
(0, 1,799) 
(0,. 1.763) 
(I), 1.740) 


ItEc^f, eoclUean nJcn of fourth order truncation error coefficients 

stability Linltss Real axis, B* e(R*, o) 

* * 

iBaginaty axis, I e(0, 1 ) 



orlor 1,KT(2)3 Al,o- 

rlthm with Third Order scaled Solution (No Additional 
Derivative Evaluations) . Acwitional 

olnen!^ ^solute value of the error in each com- 

ponent over one revolution vs. sigma. 


Seqwested Ascur-icy a 0 . WOOD-oi 
»usber Of t>,U,xtU miuatiSfis = 94 




AVj Ifiap. y I j 

^ 2 “ 


» (unscsied .solution) 

•#»9 lERP y.l. Avg IEA(^ y j_ 

3 g ■ 


0-I3M,2m, 0.14,e,5»-0, 

0.«!0Sl4^rwm a 


0.435 )45I1»01 
0. t055J6O-0t- 
0.435953H-01 
0.1061360-01 
0.1363310-01 
0. 1064480-01 
0.1065550-01 
0. 1066 5CI0-OI 
0. 1366810-^01 


0.6701800-02 
0.676512C-02 
0.6831200-02 
0.6809350-02 
0. 69686 Or-02 
0.7038330-02 
0.7107910-02 
0. 7176 5 30 -02 
0. 7243650-02 


0. 1617030-01 
0.4634450-01 
0..16S1480-01 
0..166887D-01 
0. 1686370-01 
0. 1703800-01 
0. 1720930-91 
0. 1737610-01 
0. 1761650-01 


Pequested Accuracy a 
Kunbec of Pacivi-tive ] 

ligaa A-»| jeak y_j 
1 


0. 1218470-01 
0% 1223790-01 
0. 1228990-01 
0. 1233570-01 
0.1237450-91 
0.1240580-01 
0. 1242930-01 
0..12444JD>01 
0.1248080-01 


^ . 0.10000-03 

E Pacivulve Bvalaationfe a 284 . 

(unsealed solution) 

M t... 


0.7318290-04 
0. 73-23 MD*04 
0. 7329-280-04 
0.4335470-04 
0-4343300-04. 
0.7351166-04. 
0. 735940r-04 
0. 7366520-04 
0.7375150-04 


0. 1333850L-03 
0_1342490-03. 
0..1346 160-03 
0. 135025D-03 
0.1354510-03 
”0.-l35883r-03 
0.1363244-03 
0.1367610-03 
0.1371850-01 


0.22.38980-03 

0.2250770-03 

0.2261930-03 

0-2273920-03 

0-^285970-03 

0.2298140-03 

0.2310410-03 

0.2322850-03 

0.2334970-03 


0.231014D-03 

0.2312440-03 

0.2314810-03 

0.2317320-03 

0-2319890-03 

0.2322480-03 

0.2325070-03 

0.2127330-03 

0.2324580-03 
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^ 'Pablo 3 Coofficionts for the Rungo-Kutta-Pohlborg (Fourth) Fifth 

Order Algorithm, RKP(4)S. 

I 


i 

a 

c 

c 


j 

j 

J 

0 

0 

25/216 

16/135 

1 

1/4 

0 

0 

2 

3/3 

1408/JJ565 

1408/12825 

3 

U/1J 

219V4104 

6656/56430 

4 

1 

-1/5 

“5/50 

5 

V2 

0 

2/55 


1 



k« 0 

t 

2 

3 4 

1 

1/4 




2 

3/32 

9/32 



3 

1932/2197 

-7230/2197 

7296/2197 


4 

439/216 

-8 

3680/513 

-845/4104 

5 

-8/27 

2 

•3544/2565. 

1859/4104 -11/40 


5 

Error estimate = |y-y| = | £ ((c -c.)*f.)| 

j“0 ^ J 3 

Notoj 5th order solutions are formed with 5. coefficients! 4th 
order solutions with ^ 
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Toble 4.^Charaotoristic Properties for the Fourth order, Single 
Stage, Sealed RKF(4)5 Algorithm 




0.00 

0.05 

0.10 

0.J5 

0,20 

0.25 

0.30 

0«35. 

0.40 

0.45 

0,50 

0.55 

O.OO 

0,65 

0.20 

0.75 

0.80 

0.85 

0.90 

0,95 


ITf O i 
5- 

0. 103 >24 0^02 

0J.272 130-03 
3,. 127646 0-03 
0 J00328D-03 
3. 1172440-02 
0. 1491J1D-72 
0. 172644 0-02 
7.1061950-02 
J. 1904016-02 
0, 1057600-02 
0. 1743050-02 
0. 1582270-02 
0.1396580-02 
0. 1202490-02 
0. 1003510-02 
0. M4987D-03 
3.6181600-03 
0.4176290-03 
0. 2244980-03 - 
0.6797230-04 


4t 

fi- 


* 

1 - 


(-1.A20^ ,0i. <3.274, 3.202) 


(-4.434, 0) 
(—3.622,., 0) 
(-3.241,. 0) 
(-3.016, 0) 
(-2-872,. 0) 
(-2-778,. 0) 
(-2. 722-, .0) 
(-2-697, 0) 
(-2.703^ 0) 
(-2.745, 0) 
(— 2.837, 0) 
(-1.032,. 0) 
(— 3^.786 , _ 0 ) 
(-3-4 18, 0) 
(-2.979, 0) 
(-2..801, 
(-2.738, 
(-2.774, 


0) 

3) 

0 ) 


(-3.0020 0) 


(3, 0) 

( 0 , 0 ) 

(3, 0) 

( 0 , 0 ) 

( 0 , 0 ) 

( 0 , 0 ) 

( 0 , 0 ) 

(3, 0) 

(3» 3) 

( 0 , 0 ) 

(3* -3) 

( 0 , 0 ) 

( 0 , 0 ) 

( 0 , 0 ) 

(3, 0) 

(3, 0) 

<3, 0.6093) 
(3, 0.6282) 
( 0 , 0 ) 


tT8C^(,_FucUdaaa r.oca of fifth oedoe truncation 

Stability Limits: aeal axis, a* e (R*, 0) 

. . * * ♦ 

Itt»a^in4ry akIs^ X e ^ t j 

1 2 


error coemoient 
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Table 4.b Accuracy Comparisons for the Fourth Order RKP4(5) Algo- 
rithm with Fourth order Scaled Solution (One stage) 

Average of the absolute value of the error in each com- 
ponent over ono revolution vs. sigma. 


fiegaested Accuticy « 
Noaber of Derivative 


sign a 


Avg..jBfia y I 
1 


0. 1000D-03 

Evaluations =117 (unsealed solution) 

Avg JERR y I Avg y |_Avg IBRB y | 


0.22016M0-02 0.2J634810-02 0.52830690-02 O.3862028O-O2 


0.21992400-02 
0.22J5637D-02 
0.J21 25360-02 
0.22192 320-02 
0.22254990-02 
0. 22303320-02 
0. 22339630-02 
O.2235830O-O2 
0.22360270-02 


0.16627200-02 0 , 
0. 1694677D-02 0 . 
0. 17271510-02 0 . 
0. 17601040-02 0 . 
0.17934710-02 0 . 
0.48270560^2 0 , 
O.186079fl0-02 0 . 
0.18964880-02 0 . 
0.19336200-02 0 . 


.30980840-02 

31852276-02 

.32719721>-02 

.3388099h-02 

3S117780-.02 

36347800-02 

37566500-02 

38772720-02 

39966760-02 


0.38691120-02 

OJ889998D-02 

0.39081340-02 

0-39239800-02- 

0.39397550-02 

0.39526880-02 

O.J9626200-02 

0.39693430-02 

0.39727240-02 


Reguested Accuracy. = 0.10000-05 

uabet of Oerlvative fivaiuatloas « 278 (unsealed solution) 

Si9«a Avg lEES y^, ,BRp ^ ^ ^ ^ ^ 

^ 3 4 

0.00 0.88536536-04 0.86683240-04 0.17541916-03 0.-U36550D-03 


0.83480320-04 

0.88526330-04 

0.88555150-04 

0.88601560-04 

0.68443310-04 

0.88676740-04 

0.88704300-04 

0.88720320-04 

0.83492190-04 


0.7603060D- 

0.76586186* 

0.77158040* 

0.J7706296* 

0.78199696- 

0.78711360* 

0.79244620- 

0.79764090- 

0.80229370- 


04 0.14174300-03 
>04 0.14356340-03 
04 0.14544130-03. 
04-6.14731130-03 
04 0.14910150-03 
04 0. 15088570-03 
04 0. 15266220-03 
04 0.15441410-03 
04 0.15611350-03 


0.15327480-03 
0. 15327630-03 
0.15323830-03 
0. 15317580-03 
0.15310550-03 
0.15302500-03 
0.15294760-03 
0.15289090-03 
0.15280120-03 
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Tablo 4,c Accuracy Compnrioons for tha Fifth OrdOr RKF(4)S Algo- 
rithm with Fourth Order Scaldd Solution (Ono stago) . 

Average of tho absolute value of the error in oach com- 
ponent over one revolution vs. sigma. 


Beqaested-Acaufucy 

»«abtte-a£_nj*-A 


o 0. tO0(H>-O3 
70 Evaluations 


8 


^ ■(iflscalod-soltttloft) 


Siqoa 


0.10 

0.20 

0,30 

0.40 

0,50 

0.60 

0.70 

0.80 

0.00 


ivi. teaf y I. 

1 

0.4771240-03 

0.2806)50-03 

0.2848250-03 

0.2882740-43 

0.2)33/40-03 

0.2)68300-03 

0.2994100-03 

0.3007880-03 

0.3008570-03 


Avg Ifijtp y I 
2 

0,7683450-03 

0.7751020i»03 

0.7817980>03 

0.7887860-03 

0.7960450-03 

0,8036830-03 

0.4116420-03 

0.8198160-03 

0.4280780-03 


Avq IFFP y I 
3 

0.6062570-03 
0.6329350-03 
0,6613570-03 
0.6899950-03 
0.7187510-03 
0.744956D-03 
0.7676720-03 
0. 7865270-03 
0. 80 1 5080-83- 


IBBR y I. 
4 

0.4127010-03 
0.8241490-03 
0.8364190-03 
0. 8483840-03 
0.8565500-03 
0.8664170-03 
0.8713520-03 
0.4730930-03 
4.8716570-03 


PequestM Accuctcy » 0. IOOOD-05 
Nunoer of Oaci-vatlve Evaluations = 278 


(unsealed solution) 


Siqsa 


0.04 


0,40 
0 , 20 
0.40 
0.40 
0.90 
0.60 
0.70 
0.80 
0.90 


Avj |E3P y i. 

1 

0.3/37690-04 

0.3/36160-04 

0.3/39820-04 

0.3/48330-04. 

0,3743700-04 

0.3746)90-04 

0.3749380-04 

0.3751480-04 

0.3752460-04. 

0.3748420-04 


Avg (BPE y | 
2 

0.3539990-04 

0.3053500-04 

0. 3076890-04 

0.31006/0-04 

0.3123880-04 

0.3142140-04 

0.3162950-04 

0.3186930-04 

0.3211510-04. 

0.3229940-04 


AVg LBPP y I 
3 

0, 7731490-04 

0.6234380-04 

0,6310540-04 

0.6391150-04 

0.647056D-04 

0.6540740-04 - 

0.661 0550—04 

0.6687760-04 

0.6765-780-04 

0.684391D-04 


Avg_|BRP y (_ 

8. 

0.6416440-04. 

0.6397820-04 

0,6397380-04 

0.6399110-04 

0.6400470-04 

0.6402490-04 

0.6402660-04 

0.6399880-04 

0.6397140-04- 

0.6391670-04 
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...n CoofCiolontB for tho Pivo gtago, Fifth Order, Bcalod 
RKP (4)5 Algorithm . 

Coofflclonto are Hated for oe1, with 3a < ® 04 ®0, 

for j “ 6,7,8,9,10, with a., 0. . given ln*Tablo 3 for 
J<6. ^ 


«t / 25 


228.10391 / 512000000 


3 / 1 

9,0 

3 s 

-1803163583 / 5632000000 

29 / 30 

9,3 

3. “ 

28308821 / 128000000 

39 / 40 

9,4 

3 * 

-165661179 / 352OO0O0O 

29 / iO 

9,9 

0 

1 / 2 

113121 / 781250 
-2379331 / 0593750 

9,6 

0 

9, .7 

0 

9,5. 

1 / 2 
1 / 2 

76824 / 390625 

3 

121551427 / 2400000000 

753t»96 / 4296875 

10,0 
0 » 

-8315477879 / 8800000000 

711437 / 32000000 

10,1 
3 = 

10,4 

0 

65813119 / 600000000 
-954650581 / 16SOOOOOOO 

174857013 / 352000000 

10,5 

3 » 

1 / 2 

-2604771 / 8000000 

10,6 
0 * 

1 / 2 

1245429 / 22000000 

10,7 

3 

1 / 2 

1 /'2 

10, 8„ . 
3 « 

”l / 2 


10,9 



1164791J f 300000000 
-S33t&5;63 / 1100000000 
I6t'i1543 / 75000000 
-39081107 / 206250000 



A 

0,1 

" 18644 / 

A 

0#2 

• -76165 / 10179 

A 

0,1 

• 43625 / 10179 

A 

0,4 

• -20900 / 10179 

A 

6ft 

» 1V7U2^-/ 169932 

A 

6,2 

• -721S623 / 382347 

A 

6,3 

• 7203125 / 392347 

A 

6,4 

“ -2M0000 / 382347 

A 

7,1 

“ 12064 / 2601 

A 

7,2 

» “166880 / 7803 

A 

7,3. 

“ 718000 ./ 234 09 

A 

7,4 

• -320000 / 23409 

A 

8, J 

“ -3066859635697684 / 1337382789874677 

A 

8,2 

• 274786294.95572380 / 4018148369624031 

A 

8,3 

• -58592286924705500 / 12036445108872093 

A 

8,4 

• 5444997483920000 / 12036445108872093 

A 

■ -371200 / 452907 

A 

9U- 

• 5363200 / 1358721- 

A 

9,3 

« -25120000 / 4076163 

A 

9f4 

• 12800000 / 4076163 

A 

10,1 

» -5795511385216816 / 133M82789874677 

A 

10,2 

• 88186674341465120 / 4012148369624031 

A 

10,3 

■ -387934866359732000 / 12036445108872093 

A 

10,4 

• 176347228586080000 / 12036445108872093 
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Table s.b CharaotorlstiQ Proportion for tho Fifth Order. Pivo 
Stage, fiealod RKP(4)S Algorithm. 


Slgaq 

0.00 

o.os 
0.10 
o.i<i 
0.20 
0.2S 
O.JO 
0.39 
0.40 
0.4 S 
O.SO 
0.55 
0.60 
0.65 
0..70 
0..75 
O.flO 
0.05 
0.40 
0..95 


iJlSC I 
6 

0. 31 )5/75to02 

o.aiissoio-oi 
3. 13310(t(i0«02 
■).32J0340U-12 
3. 2'3 305140- )2 
0. 13320430-02 
0.34304600-32 
0.32 i34640L>-02 
3. 29 3 1 1420-02 
0.27147750-02 
0.27178590-02 
0.27361370-02 
0.27195630-02 
0. 26323280-02 
0. 22112260-02 
3. 21 466680— 02 
0. 25175490-02 
0. 31/92050-02 
0. 3 4133830-02 
0.25209160-02 


(-1.678, 0) (2.046,3.607) 


(- 1 . 202 ,_ 0 ) 
(-3. 128, 0) 
(-3.113, 0) 
(—3-131, 0) 
(-3.189„ 0) 
(-3. 340, 0) 
(-3.146, 0) 
(-3.049, 0) 
(-2.899, 0) 
(-2.809, 0) 
(.-2.730, 0) 
(-2.663, 0) 
(-2.602* 0) 
(-2.944., 0) 
(-2,49$, 0) 
(-2.461, 0) 
(-2.466, 0) 
(-2.951, 0) 
(-3.012, 0) 


( 0 . 0 ) 

( 0 . 0 ) 

0 , 0 ) 

10 , 0 ) 

( 0 , 0 ) 

( 0 , 0 ) 

(t.sot, 2 


( 0 , 

(0» 

(0# 

( 0 , 

( 0 , 

( 0 , 

(0. 

(0 


( 0 , 0 ) 

( 0 » 0 ) 


2-683) 

2*410) 

2.219) 

2 . 080 ) 

1*965) 

1.846) 

1.638) 

0 ) 

0 ) 


.981) 


•TTC^I, RucUdoaa nota of atuh otdor ttuncatloa-^rtot coefficients 
stability Liaits: 6cal .tala, 0) 

1 2 
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■f«Wo G.o Aoournoy Compnrinonfl tho Fifth Ordor RKP(4)9 Al«o- 
vlUm with Fifth ordop floalnd flolution (Five Gtagos) . 

Avoro.jo of tho abooluto valuo of tho orror in oaoh oom- 
ponont ovar onn revolution vn. oiymn. 


Pt>/|iionto(l Auuui i«y n 
•JiJobAr of OoiiviuvL' 


0. IOOOD-0,1 

Sva luett lotui o 


117 


(unticalad uolutlon) 


Pl<l«a AV't l"l'F. y | 
1 

^•00 0. 2MJ 1l)0a>0j 

0. 10 0.202'iM3-03 

0.20 0.2906605-03 

0,30 0.2802513-03 

0.«0 0,291J07)-03 

0.50 0.29)K*2J-0J 

0.60 0.2927383-03 

0.70 0.28)7131-03 

0.80 0.26)3570-03 

0.90 0. 2818830-03 


8tiqaost(>8 Auuncy « 
Nunbet; of Oiclvjttivc 


8iqna AVq ]£Rt y I 
1 

0.00 0. 37timM)-04 


AVI I ERR y I 
2 

0.9016930-03 

0, 7801640-03 
0.7979970-03 
0.8084920-63 
0.8130820-03 
0.3174660-03 
0.a28221l'-03 
0.8495220-03 
0.8798550-0.1 
0.9096 040-03 


0. 10000-05 

E9a4.uatioi,a 


A»9 IE38 y | 
2 

0.3540930-04 


AV(J lERfi y ( 
3 

0.9877650-01 

0,6342610-03 
0,6879770-03 
0.729663P-03 
0. 7561380-03 
0.7751540-03 
0.7982580-03 
0.8311530-03 
0.9116520-03 
0,9843970-03 


A»9 IPRp y I 

3 

0.7740810-04 


Avq I EBP y |_ 
4 

0.8170480-03 

0.8234080-63 
0.84.35500-03 
0. 8520170-03 
0.8476046-03 
0.8359580^03 
0. 8252910-03 
0.8220620-03 
0,8.116870-03 
0.8407440-03 


solution) 


Avq IfieS y I 

4 

0-6422690-04 


278 (unsealed 


0. 10 
0. 20 
0.30 
0.40 
0.50 
0.60 
0.70 
0.80 
9.90 


0i374693i)-04 
0. 175)8)^-04 
0.37)5723-04 
0. .1746853-04 
0.37)36)0-04 
0.373)1'ii)-04 


0.307795D-04 
0.312677D-04 
0.3179710-04 
0.3233650-04 
0,3289300-04 
0.3341506-04 
0 


0.63-24040-04 

0.6478210-04 

0.6626990-04 

0.6778870-04 

0.6927970-04 

0.7084270-04 

0,7254550-04 

0.7437860-04- 

0,7610000-04 


0.6422090-04 
0. 6443200-04 
0.6448590-04 
0.6440720-04 
0.6437540-04 
0.6418760-04 
0.6416790*04 
0.6422630-04 
0. 6423680-04 


0.3n655j-04« 0.3433^51)- 
0. 373 J 593-04 0*34fl043D- 
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Table (5.a ?roo P«ran.eter eolootion ana Charaotorlotle Proportion 

for tho Fifth Ordor, Two fita^o, Sralod RKF4(0) Alaorlthm. 




rec 


o*rto 

o.sft 

0,10 

OvW 

0^20 

OviS 

0.30 

0,15 

o^to 

0.49 

O.M 

0.55 

0.60 

0.65 

0,70 

0.75 

O.flo 

0.65 

6.90 

0.99 


0.J510 
0„2500 
0.4500 
0.4166* 
0.25P0 
0.3750 
0.1500 
0« 3 1 25 
0.2727* 
0.5200. 
0.2600 
0.3300 
0.4500 
0..4600 
0.4500 
0.4300 
0. 42 00 
0.4100 
0. 3750 


0.650 
0.800 
0.9 10 
t.lOO 
1.000 
t.OOO 
1.010 
1.000 
1.000 
6.4 70 
0.630 
0.720 

oaoo 

0.590 

0.590 

0.590 

O.SSO 

0.500 

O.J)-50 


1.3.1560-02 ( 1.678,0) (2.045,1.607) 


1.25460-06 

1.45130-05 

0. 1 1620-03 
3. 176 70-04 

1. 16420-02 
9,23470-03 
1.49770-03 
1.85750-03 
1. 15000-02 
0.32950-03 
0.32210-02 
0.38120-02 
1.17750-02 
0.. 14300-02 
0. 36510-02 
1. 37700-02 
0,36110-02 
0.36090-02 
0. 19100-01 


tepeatina doclaatu 


(- 11 . 12 '»* 0 ) 
( -7.031,0) 
( -4.110,0) 
( -5.311,0) 
( -2.511,0) 
( -4.140,0) 
( -1,514,0) 
( -3.416,0) 
( - 1 . 21 . 1 , 0 ) 
(. -3.402,0) 

( -2.765,0) 
( -2,700,0) 
t -1.210,0) 

( -1.997,0) 

( -5.202,0) 

( -2.701,0) 

( -2.789,0) 

{ -2. 270,0) 

{ -1.945,0) 


( 0 , 0 ) 

(1,3.604) 

( 0 , 0 ) 

( 0 , 0 ) 

(9,2. 520) 
(3-.2i3,4. 104.) 
(2.814,3.492) 
(2.516,1.581) 
(2.160,3.407) 

(0,.1.1fl1) 

(1.988,2.991) 

(1,870,2,971) 

(2.108,3.616) 

( 0 , 0 ) 

( 0 , 0 ) 

(1.605,1. 120) 
(1.605,3.232) 
(0,1.142) 
( 1 . 260 , 2 , 208 ) 


Ptaa -Parameter's 


a , 

6 


Chatadtectstia lUrattEjctLaas 


StmiHy LlnUss Beil axlu,.. 8* e (R*, O) 

* <1 III 

lacy ax-ia, i e (i , i ) 
1 2 
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TJjblo 6.b Coefflclortts for tho Two 8togO| Plfth-Ordor*. 8oalod~ 
RJ<s^4)5 AlgorlthJh for Valuos of tho Scaling Baramotor, 


M jl » a /a 

/ o 

j#lL 


k®0, 1,,.,, j-i (f COB Table 


CO 

c< 1) 

C( 2) 
C( 3) 
CC 0) 
C< 5) 
C( 6) 
C( 7) 

BO (6) 

B(6,1) 

B(6*2) 

B(6,l} 

B(6,4) 

B-U*5) 

BO (7) 

B(7,n 

B(1»,2) 

B(7,3) 

B(7,t) 

B(7,5) 

B(7»6) 


0. 2000000000000000D«00 

0. 41 6666666606606 70^00 
0. I'MOOOOOeoOOOOOOD^Ol. 

0. 12496352330962360400 
0. 00000000000000000400 
-0. 174.12801118442410-01 
-0. 1 7266412915603870-03 
0. 82891774891770620-04 
0. 4 130 724386 72433 10-02 
0. 62687602569482260400 
0. 26 U062995816362iUOO 

3. 32545914297783570400 
0.1 32552646550.1 1150400 
0. 0'3000000000000000400 
-0..3 38839651832299.10-01 
0.23538942491949290-01 
0. 00000000000000000400 

-0.24252117427572310400 
-0. 3 J82978320350067D400 
0-1 1871822499691400400 
0.57266104666808820-02 
-0.90433530513109870-02 
0.4 7066 19454105a 000-01 
0. 11983513293573920401 


A(6) 

M7) 

60 

C( 1) 
c< 2) 
C( 3J 
C( 4) 

SI 

C( 6) 
CC7) 

B0<«) 

B{6,2) 

B<6,3) 

B(6,4) 

B(B^5) 

BO (7) 

B(7, 1) 

9(7,2) 

B(7,3) 

B<7,4) 

9(7, S) 

B(7,6) 


4 ( 6 ) 

4(7) 

CO 

C( 1) 
C( 2) 
C( 3) 
C( 4) 
C( 5) 
C< 6) 
C( 7) 

BO(JS) 

B(6,t) 

B(6,2) 

B(6„3) 

B(6,4) 

B(6,S) 

BO (7) 

9(7, 1) 

8(7,2) 

8(7, 3) 

8(7,4) 

8(7,5) 

8(7,6) 


• O.43OOOOOOOOOOOOOO0»OO 

■ O.J1280000000000000*00 
» 0. UOOOOOOOOOOOOOODfOl 

• 0.9O98835070834918O-O1- 

• 0. OOOOOOOOOOOOOOOO&^OO 

• 0* 19116897666368260*01 

' ”j*-l23i2042ll 1949700-02 
' 7233353760*03 

' 0.9218935^669356760-01 

5*?®!®?(3273626358n*00 
-0.11 840675664203070*0 1 

0.19974879930870760*00 
8* 19489306819995050*00 
0. OOOOOOOOOOOOOOOOOtOO 
-0..1 J43230 583627 1650*00 
0. 8208 1190854058330-0 1 
0, 03000000000000000*00 

0.24002487612025570*00 
0. 11420294634216480*01 
26922744501220540*00 
®**1 212229 7761278770*00 
*0* 73168790140208820-01 
-0. 7744341 8043669970-01 
-0.28442766395860750*00 

* 63000000000000000*00 

* ®- 93000000000000000*00 
“ 9.72000030003000000*00 

“ 9* 99560344717122290-01 

* 9*09000000000000000*00 
» -0.64892741212537400*00 

* M??21798325569660-01 
“ “9.71136124276835130-02 
» -9.24774011808635590-92 
I ?*512S6598353146880*00 

* 0..9)75702931597602D*.00 

I ?*’J?2d26l 139154680*00 
I ^ ^2469110602688580*00 

* ®-93990000000000000*00 
. ”2*'?«*i?-^23 185224540*00 

^♦17740851443381280*00 
» 9.0)000000000000000*00 

“ Mt!!7!f2192468717SO*00 

* ,?• l'i?^177657728288D*00 

“ 27257766742750820-02 

r ‘2* 595250-01 

“ 29019378869206320-01 

» ?* |!294261579371260*00 

9. 222 162 14626640500*00 


4(6) 

M7) 

■ CO- 
C( 1) 
C( 2) 
C( 3) 
C( 4) 
C( 5) 
€( 6 ) 
C( 7)_ 


80(6) 

B(6, 1) 

B(6,2) 

B(6,3) 

B(6,4) 

B(6,5) 

80(7) 

B(7,1) 

B(7,2) 

M7,3) 

B(7,4) 

B(7,5) 

B(7,6) 


■ 9» OM00000000000900*00 

® 9.43009000000000000*00 

• 9.53000000000000000*00 

• ®*2!I^^®7‘*980993280«*01. 
® 9.09000000000000000*00 

■ -9.819991115S151296D*01 

■ 9.11909993304736690*00 

• 7796500702063710-01 
! "M!I??®94S1958462D*00 

• ®*®791 1256913501890*01 
-M_ 3. 3713t934d527065.70*01 

• 9.13215873920594960*00 
- 9.23647154938848580*00 

• 9. ojooooo6ooooooooo*oo 

• ?*’i®587837d49l868D-01 
« -9, 18489072379354070-01 

• 9.09000000000000000*00. 

• “S* 2799871 0662899930-02 
■ 2* !I?§7?22826 333870*00 
I “2*|??I®737636920460*00 

71 13607481 1662960-01 

• “9.30026854561651860-01 

• 9.83766740264807390-02 
“ 9* 85496024395663850-01 
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table Aeburaoy Comparisons for the Fifth Crder RKF(4)5 Algo- 
rithms With Fifth Order Scaled SoluUort (Tw6 Stages) 

Average of the absolute value of the error in-^each oo«n> 
poAent over one revolution— vs . sigma. 


ffeqoested Acoucicy = 0.10000*0-3 

iluabec of Deciv4ti«e-F»aitutiOfts « 117 tunscaled 

Sigsa Avi; |38B y |_.A»g y I Avg. JfiSR y I 

1 2 ^3 

rt-OO 0.2800020*03 0.9030270-03 0.987S20D-03 


solution) 

Avg |8FR y J 
4 

0.8167670-03 


0.10 

0.20 

0.30 

0.40 

O.SO 

0.60 

0-70 

O.S0 

0.JL0 


0.2768340-03 

O.277*240r03 

0.278426D-03-- 

0.279 >260-03 

0.2745230-03 

0.2845380-03 

0..271525D-03 

0.-2778120-03 

0-2764150-03 


0-7760320* 
0.7907130- 
0.8054480- 
0.8212740- 
0.33-29820- 
0.8547150- 
0.8609340- 
0.8771070- 
0. 8846570- 


03 0.6212980-03 

•03 0.6580050-03 
03 0.69565-70-03 

03 0.7400680-03 

03 0.-7576860-03 
03 0.8365420-03 
03 0.80815)0-03 

03 0.8558320-03 

03 — 0. 8962210-03. 


0.8158670-03 
0-8248950-03 
0. 8312730-03 
0. 8378940-03 
0.8309540-03 
0.8570090-03 
0-8202220-03- 
0. 8288250-03 
0. 77995*0.-03 


Requested* Ascur.icy « O.l6o0o^O5 

Buabet »f Daciyitlve Evaluations = 27.8 (unsealed solution) _ 

-Eigna- 

Avj (EBF y.. 1 
t 

Avg (Pftp y ( 
2 

Avg lEFft y (_ 
3 

Af9 t l_ 

4 

0.00— 

0. 3737690-04- 

0.3539990-04. 

0.7731490-04— 

0.64 16440-04. 

0,10 
0.20 
0.30 
0.4 0 
0.30 
0.J6O 

0.70 

0.80 

0.90 

0. 3736630-04- 
0.373 7130-04 
0.373 >540*04 
0J74171O-04- 
0.3733720-04 
0.3743660-04 
0.3736900-04 
0. 3-747810-04- 
0.371453 0-04 

0.3078780-04 
0.3130390-04 
0.3174960-04 
0.3229660-04 
0.3272780-04 
0. 3333530-04 
0.3379790-04 
-0.34-2842D-04 
0.3431610-04 

0-6305490-04 

0.6460140*04 

0.6594550-04 

0.6767080-04 

0.6912960-04 

0.7109360-04 

0.7230900-04 

0-737902D-04 

0.7547140-04 

0.6408260-04 

0. 6420190-04 

0.6428430-04 

0.6436620*-04 

0.6417080-04 

0.6432790-04 

0.5406600-04 

0.6394770-04 - 

0.6.165250-04 
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